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One possible low energy process due to beyond the Standard Model (BSM) physics is the neutron- 
antineutron transition, where baryon number changes by two units. In addition to providing a 
source of baryon number violation in the early universe, interactions of this kind are natural in 
grand unified theories (GUTs) with Majorana neutrinos that violate lepton number. Bounds on 
these oscillations can greatly restrict a variety of GUTs, while a non-zero signal would be a 
"smoking gun" for new physics; however, to make a reliable prediction, the six-quark nucleon- 
r- 1 antinucleon matrix elements must first be calculated non-perturbatively via lattice QCD. We re- 

view the current understanding of this quantity, describe the lattice formalism, and present pre- 
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liminary results from 32 3 x 256 clover- Wilson lattices with a pion mass of 390 MeV. 
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1. Introduction 

One unanswered mystery of the universe is the process that led to the abundance of observed 
baryons as compared to their antibaryon counterparts. The source of this baryon number violation, 
which is expected to come from beyond the Standard Model (BSM) physics, can be realized in 
low-energy processes such as proton decay (if baryon number is violated by 1 unit) or transitions 
between neutrons and antineutrons (if baryon number is violated by 2 units). The latter case, 
often referred to as neutron-antineutron oscillation (akin to neutral meson mixing), proves to be an 
intriguing scenario when considering the usual sphaleron picture of baryogenisis (which violates 
baryon number, B, and lepton number, L, but conserves B — L) coupled with Majorana neutrinos [1] 
(whose transition between neutrinos and antineutrinos leads to AL = 2). Additionally, neutron- 
antineutron oscillations do not suffer from kinematic suppressions that can restrict proton decay if 
there is little overlap with the initial state proton and final state electron or muon [2]. To that end, 
neutron-antineutron oscillations have been explored experimentally with intriguing prospects for 
^(1000) improvements in upcoming experimental efforts [3]. 

Any discussion of neutron-antineutron oscillations starts with assuming the existence of some 
BSM process that leads to a AB = 2 operator in the low-energy effective field theory. This operator 
will lead to off-diagonal elements of the Hamiltonian of the neutron-antineutron system 

ySm En) \ 8m E-VJ 

where V is the potential difference between the neutron and the antineutron (a magnetic field can 
lead to a non-zero V since the magnetic moments have opposite signs) and V = in a free system. 
Upon solving the Schrodinger equation for the system, one finds the transition probability between 
neutrons and antineutrons given by 



(1.2) 



While this equation is true for a given V, it is standard to define the period of free neutron oscilla- 
tions due to the BSM physics as 

*m=^-- (1.3) 
dm 

The value for x m greatly depends on which BSM scenario is being explored. It has been estimated 
that a bound of x nTt > 10 10 — 10 11 seconds is sufficient to rule out many of the current models [4]. 
For example, TeV-scale seesaw mechanisms for neutrino masses in SU(2)l x SU(2)r x SU(4) c 
are expected to be ruled out at % m > 10 10 — 10 11 seconds [5]; and 50(10) seesaw mechanisms 
with adequate baryogenisis, at x m > 10 9 — 10 12 seconds [6]. Current experimental limits can even 
restrict extra-dimensional models with new particles with masses below a TeV [2,7]. It should be 
emphasized that these are order of magnitude estimates with the QCD input coming from naive 
dimensional analysis. Future estimates will require rigorous and precise lattice calculations to keep 
pace with experimental precision. 

The detection mechanism for these transitions is the cold annihilation of a newly formed an- 
tineutron with a nearby neutron (for more details, see W. M. Snow's plenary at PXPS 2012 [3]). 
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The primary channel for this cold annihilation is nn — >■ 5k, and this unique signature allows for 
experimental signals with little or no background. Generally, there are two sets of experimental 
searches. The first, which comes for free with large proton decay detectors such as Super-K, is 
based on neutron-antineutron annihilation within nuclei. Naively, one might expect this to occur 
quite frequently, as the number of nuclei far exceeds the expected bound x m > 10 11 ; however, the 
oscillation period within nuclei is highly suppressed, with a magnitude of roughly [8] 

% uc / = (3xl0 22 )^. (1.4) 
sec 

As a result, the bound is suppressed compared to the free expectation, and one must rely on model 
estimations and extrapolations to extract it. To date, the most stringent bound from experiments of 
this kind, % m > 3.5 x 10 s seconds, comes from Super-K (2011) [9]. 

The second type of experiment explores the annihilation of free, cold neutrons with a target 
after a significant time of flight. This type of experiment is free of the model-dependent estimations 
required for annihilations within nuclei and allows for greater control of systematics. To date, the 
most stringent bound comes from the ILL experiment (1993): T nH > 0.86 x 10 s seconds [10]. A 
factor of ^(1000) increase is estimated for future experiments of this kind, but the bounds to 
rule out various BSM theories could be altered significantly depending on QCD enhancement or 
suppression of the neutron-antineutron matrix elements. 

2. Oscillations and matrix elements 

The observed value of the mixing arises from three inputs 

— = 8m = c B sm(I*bsm, 1*w)cqcd(Vw ,AQCD)(n\<^\n) , (2.1) 

where cbsm is the running of the BSM theory to the weak interaction scale, cqcd is the QCD 
running from the weak to the nuclear scale, and (n|^?|n) is the non-perturbative matrix element 
mixing the neutron and antineutron states. The one-loop perturbative QCD running, cqcd, is known 
[7, 11], and cbsm has been calculated for multiple theories [2,5-7]. The operator & contains two 
up quarks and four down quarks and is composed of three pairs of quarks from the possible forms 

u T Cu , u T Cd , d T Cd , (2.2) 

where C is the charge conjugation matrix. Additionally, these terms always come in chiral pairs, 

uJ^Cdi , u\CdR , (2.3) 

since the mixed chirality terms are zero. Lastly, these operators are invariant under color symmetry, 
SU (3) c , which leads to two color tensors 

^ijklmn ^mik^njl ^nik^mjl ^mjk^nil ^njk^mih ^tjklmn ^mij^nkl ^nij^mkli (2-4) 

where i,j,k,l,m,n are color indices. These three conditions lead to three types of operators [12]: 

®X\XiX3 = ( u ixi Cu jX^(dkx 2 C dlX2)(dLx3 C dnx 3 )r S ijklmni 
®X\XiXl = ( u fxi Cd jXi)( U kX2 Cdl X2)( d liX3 Cdn X3)^'ijklmn^ 

^XiX2X3 = ( U fxi Cd JXi)( U kX2 Cdl X2)( d Lx3 Cdn X3)^ijklmn^ ( 2 - 5 ) 
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where %i = L,R. At first glance, there would appear to be 24 independent operators, but there are 
several additional symmetries. The first set of symmetries due to the flavor structure is 

*U = ^4 3 = ^X3> ^.6) 

which reduces the set to 18 independent operators. An additional symmetry that emerges from 
antisymmetrizing pairs of epsilon tensors over four indices leads to (with a = L,R and p = L,R) 
[13] 

^L p -Ka P =^lap, (2-7) 

which reduces the set to 14 operators. In addition to enforcing SU (3) c , it is also expected that the 
operators should be invariant under SU (2)l<8> U(1)y- This gauge symmetry and the symmetries in 
Eq. (2.6) leave only six operators: 

<? 4 = 20 3 LRR , & 5 = 40[ LR , & 6 = 4(0hx-0hx)- (2-8) 

We will present results for these operators; however, including the symmetry in Eq. (2.7) leads to 
the conditions, 

^ 6 = -3^5, &2-&l= 3=^3, (2.9) 

which reduces the number of independent operators to four. We will use these last two equations 
for a consistency check of our calculation. 

3. Lattice formalism and contraction details 

The mechanism to extract the neutron-antineutron matrix elements follows the common prac- 
tice of taking ratios of three-point to two-point correlation functions. In particular, the three corre- 
lation functions of interest and their large Euclidean time behavior are given by 

C m (t) = (N(t)N(0)) -+ \(N\n)\ 2 e- m »>, C m (t) = (N(t)N(O)) -+ \ (N\n) |V m » f , 

C mN (h,t 2 ) = (A^(? 2 )^(0)F(-? 1 ))^(n|]V)(^|«)e- m »( ri+ ' 2 )(«|^| ? i). (3.1) 

The desired quantity of interest, (n\ff\n), is the long time asymptote of a combination of these 
correlation functions, 



CNN{h)Cfjf}{t2)Cjjfj{ti -K2) 
CjWvl)CNN\h)CNN\t\ +t2) 



{n\ff\n). (3.2) 



The six-quark neutron-antineutron three-point correlation function has several key advantages 
over typical bi-linear or four-quark nucleon operators. First, if the starting point for the propagator 
is at the operator insertion (as shown in Fig. 1), only one propagator is needed per measurement, 
whereas the typical nucleon three -point function requires two propagators, one that starts from the 
source and one that starts from the operator. Second, because the propagator starts at the opera- 
tor, one can acquire all the source-operator separations (given by t\ in Fig. 1) and operator-sink 
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Figure 1: Comparison of neutron-antineutron three -point contractions (left) to typical bilinear three-point 
contractions (right). One propagator is required for a measurements at all (t\,t2) for the left diagram and 
two propagators are required for one measurement at a single t\ -value on the right diagram. 

separations (given by ?2 in Fig. 1), which allows for a two-dimensional analysis to quantify the 
excited state effect. Alternatively, typical three point functions require far more computational re- 
sources to quantify excited state effects. Lastly, the neutron-antineutron matrix element contains no 
disconnected or quark loop contractions, which removes the need for costly all-to-all propagators. 

One possible disadvantage is that multiplying six propagators together (as done for the neutron- 
antineutron correlator) could increase the signal-to-noise degradation as compared to bilinear or 
four-quark matrix elements; however, we find a reasonably good signal-to-noise ratio, as shown in 
Fig. 2. 

4. Lattice Details 

The lattice calculations were performed with Chroma [14] using the 32 3 x 256 anisotropic 
clover- Wilson lattices defined in Ref. [15] with a pion mass of 390 MeV. The temporal and spatial 
lattice spacings are roughly 0.035 and 0.123 fm, respectively, and the total spatial extent is roughly 
4 fm {m n L ~ 7.8). For this preliminary calculation, we use a total of 159 configurations, each 
separated by 4 trajectories, to calculate 7268 propagators with Gaussian-smeared sources. Con- 
tractions of these propagators lead to the same number of measurements at all source-operator and 
operator-sink separations. 

5. Preliminary results 

The desired matrix elements, (n\ &i\n), can be extracted from the long Euclidean time behavior 
of Eq. (3.2). For each ratio there are two time inputs, the source-operator separation (?i) and the 
operator-sink separation fe). In Fig. 2, for the £?\ operator is plotted against ?2 for six different 
values of t\. Two features stand out from these plots. First, there is a significant range of time 
slices where a signal can be extracted and the signal-to-noise degradation is not overly restrictive. 
Second, it is evident that there is significant excited state dependence as t\ is varied (for example, 
the plateaux extracted for t\ = 10 and t\ = 30 are significantly different). For this reason, it is very 
important to use all information available to explore the full behavior of & as a function of both t\ 
and ?2- 

In Fig. 3, the 2D plot of Si and \3t\ are plotted against t\ and ?2- Again, it is clear that there 
is a significant amount of non-trivial behavior due to excited states. To that end, a 2D correlated 
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h = 5 ti = 10 ti = 15 




h = 20 t x = 25 ti = 30 




Figure 2: Plots of ^ vs. ti for six values of t\ — 5, 10, 15,20,25,30. The large t\ and ti behavior of ^ 
should approach the neutron-antineutron matrix elements of interest. For t\ > 30 (t\ > 1 fm), the plateau 
does not change appreciably, but signal-to-noise decreases. 



fit has been performed over the time slices 10 < ?2 < 25 and 30 < t\ < 40. For this preliminary 
calculation, systematic errors are estimated by adjusting 2D fit window ±1 on all sides. 

The bare (unrenormalized) results for (n\&i\n) are shown in Table 1. Eq. (2.9) is satisfied 
exactly, configuration by configuration, but only stochastically in Table 1 due to the bootstrapping 
in the analysis. The corresponding values calculated from the MIT bag model are also displayed for 
comparison. The magnitude of each operator as computed on the lattice is below that derived using 
the MIT bag model; however, it should be emphasized that the lattice results are very preliminary 
and still require renormalization factors. 




Figure 3: Two-dimensional plot of (left) and \&\ (right) as a function of t\ and t% ranging from time 
slices to 60. The lighter colors represent larger values of ffl. The two-dimensional plateau is achieved for 
10 < h < 25 and 30 < t { < 40. 
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Operator 


Lattice Calculation (1(T 5 GeV) 


MIT Bag Model Calculation (10~ 5 GeV) 


(n\&>\\n) 


1.57 ±0.85±{j:3o 


-6.56 


(n\^ 2 \n) 


-0.20 ±0.14+^ 


1.64 


(n\^ 3 \n) 


-0.24 ±0.26±#$ 


2.73 


(n\^4\n) 


-0.02 ±0.39±g:« 


-6.36 


(n\^ 5 \n) 


0.34 ±0.82+°;^ 


9.64 


(n\&> 6 \n) 


-2.07 ±1.10+^ 


-28.92 



Table 1: Table of results for neutron-antineutron matrix elements from the bare, non-renormalized lattice 
calculation and MIT bag model calculations [12]. 

6. Systematic effects 

The primary systematic uncertainty in comparing the results of Table 1 to experiment is the 
unphysically large pion mass used. While it is not clear that an IR quantity such as the pion mass 
should dramatically effect the short distance six-quark vertex, it is a distinct possibility given that 
contractions of this system are similar to those for low-energy AW scattering, where physical quark 
masses are expected to lead to a dramatic increase in the scattering length [16]. 

The second source of systematic uncertainty is the lattice cutoff (i.e., discretization) and match- 
ing the lattice regularization to the usual MS scheme used in the perturbative running [7, 1 1] . Gener- 
ically, operators of interest might mix with lower dimensional operators with the same symmetries, 
leading to diverging 1/a corrections (where a is the lattice spacing). However, this is not an is- 
sue for these operators since the lowest dimension operator that can lead to a AB = 2 interaction 
requires six quarks. Regardless, there are expected to be Gia) corrections and renormalization 
coefficients that should be quantified. 

The third systematic which was clear from Fig. 2 and Fig. 3 is excited state contamination. 
The calculation of the six-quark neutron-antineutron correlator gives us a unique view of these 
contaminations as a 2D function in t\ and t 2 , which is difficult to come by for any other nucleon 
three -point function. For this reason, we should be able to accurately quantify these contamina- 
tions. Finally, finite volume effects should be quantified as well, though their impact is expected to 
be insignificant given the m K L ~ 7.8 lattice used. 

7. Future prospects 

We are in the process of taking several steps to improve upon this very preliminary work: we 
are extending the calculation presented here as well as repeating it for a lighter, 240 MeV pion 
mass at the same volume and for the current pion mass with a smaller, 2.5 fm spatial extent with 
a larger ensemble. We are exploring perturbative and non-perturbative lattice renormalization to 
properly match onto the perturbative QCD running previously calculated. We are also refining our 
analysis procedures to better quantify excited state effects. 

Within the next year or two, we hope to carry out this calculation both with physical pion 
masses and with a chiral fermion discretization (domain-wall fermions). Both calculations are 
numerically expensive, but within reach of the LLNL 20 PetaFlops Sequoia BG/Q. 
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